Short-Range Ising Spin Glass: Multifractal Properties. 
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The multifractal properties of the Edwards- Anderson or- 
der parameter of the short-range Ising spin glass model on 
d = 3 diamond hierarchical lattices is studied via an exact re- 
cursion procedure. The profiles of the local order parameter 
are calculated and analysed within a range of temperatures 
close to the critical point with four symmetric distributions 
of the coupling constants (Gaussian, Bimodal, Uniform and 
Exponential). Unlike the pure case, the multifractal analy- 
sis of these profiles reveals that a large spectrum of the a- 
Holder exponent is required to describe the singularities of 
the measure defined by the normalized local order parameter, 
at and below the critical point. Minor changes in these spec- 
tra are observed for distinct initial distributions of coupling 
constants, suggesting an universal spectra behavior. For tem- 
peratures slightly above T c , a dramatic change in the F(a) 
function is found, signalizing the transition. 

05.50.+q, 75.10.Nr and 64.60.A 



I. INTRODUCTION 

The understanding of the nature of the spin-glass 
(SG) condensed phase in real systems has been chal- 
lenging many authors |IJ], since the scenario emerging 
from Parisi's mean-field solution of the Sherrington- 
Kirkpatrick (SK) || model came out. Some raised con- 
clusions, like the structure of the free-energy barriers 



corresponding to many distinct phases below T c (pure 
states) arranged in an ultrametric structure and the 
existence of a critical ordering field for the condensed 
phase ||, generated controversies that remain not sat- 
isfactorily elucidated. In particular, the domain-wall 
phenomenological scaling approach (droplet model) dis- 
miss the SK model as appropriated for the description 
of short-range Ising spin glasses in low dimensions and 
does not share the same conclusions On the other 

hand, recent works based on numerical simulations pre- 
sented results indicating that short-range models should 
exhibit the same qualitative features appearing in the SK 
model H . It is worth to mention that many efforts have 
been devoted to investigate exactly-solvable short-range 
SG models as an attempt to describe real spin glasses, 
where certain aspects of the system, e.g., the correlation 
length, the sensibility to the boundary conditions, and 
finite-size effects should present a very distinct behavior 
from those of infinite-range models. For the SG model 
on the pathological Bethe lattice with finite connectiv- 
ity, the controversy about the nature of the condensed 
phase still persists. For instance, it was found that a 
replica-symmetric solution is stable for zero field when 
open (uncorrelated) boundary conditions are considered 
|)j , while the breaking of replica symmetry is required to 
obtain a stable solution below T c when closed (correlated) 
boundary conditions are imposed to the system ||lo[| . An- 
other line of approach in the study of short-range SG 
behavior was developed after the work of Southern and 
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Young who succeeded to obtain, by using the Migdal- 
Kadanoff renormalization group (MKRG) scheme, phase 
diagrams showing the presence of a SG phase in three 
dimensions (d = 3), but not for d — 2, indicating that 
the lower critical dimension di should lie in this interval. 
This latter approach, which can be viewed as an approx- 
imation for real systems, was applied to investigate the 
exponents required to describe the transition to the con- 
densed phase. On the other hand, based on the scaling 
theory, it was found that SG obeying symmetric distri- 
butions is characterized by four independent exponents, 
the thermal and the chaotic ones at T = and T = T c 
p2[ . Within this approach the chaos exponents govern 
the sign-changing of the effective coupling of two spins, 
a distance L apart, with the temperature. The interpre- 
tation for this phenomenon of critical chaos was done in 
the framework of the droplet theory [fl3f . 

Another point of view of the MKRG was explored since 
it was proved that its renormalization-group equations 
are exact for the Ising model on a family of diamond 
hierarchical lattices (DHL) [O. The study of spin sys- 
tems on such exotic lattices, whose coordination number 
varies from two to infinity, acquires relevance because 
exact solutions can be obtained and are well-controlled. 
Motivated by this, we generalized the method proposed 
by Morgado et al. jl5| (MCC) to investigate the struc- 
ture of the local Edwards- Anderson (EA) order parame- 
ter of the Ising SG model on these lattices. Through this 
method, we are able to calculate the local magnetization 
at each site by an exact recursion procedure and further- 
more, the local EA order parameter associated to each 
site. The MCC method was firstly applied to the pure 
ferromagnetic Ising model on a two-connected DHL, re- 
vealing the multifractal structure of the local magnetiza- 
tion at the critical point jig ]. Moreover, they showed that 
there is an infinite set of (3 critical exponents for the local 
magnetization which are linearly related to the a-H61der 
exponent of the F(a) multifractal spectrum. This work 
was further generalized for a p-connected DHL with q in- 
termediated sites within each connection jl6| . The MCC 
method was also applied to investigate the multifractal 
and critical properties of the q-state ferromagnetic Potts 
model on the p-connected DHL [jlTj . 

In the present paper we extend the preliminary results 
obtained by Coutinho et al., jl8| where the F(a)-spectra 
of the SG Ising model with an initial Gaussian distri- 
bution was studied for lattices up to twelve generations. 
Now we consider four distinct initial distributions of cou- 
plings (Bimodal, Gaussian, Uniform and Exponential), 
with lattices up to sixteen generations, investigating the 
temperature dependence of the -F(a)-spectru:m around 
T c . In section II, the model Hamiltonian is presented 
and the corresponding critical temperatures associated 
with each distribution are obtained. In section III, wc 
give full details of the generalization of the MCC method 
for the SG case, as well as about the numerical procedure 
used to obtain the EA order parameter profiles for the 
four distinct initial distributions of couplings. In section 



IV, the multifractal properties of the EA order param- 
eter are obtained and the corresponding F(a)-functions 
calculated by distinct initial distributions are compared. 
Furthermore, we investigate the temperature behavior of 
the boundaries of the F"(a!)-functions, above and below 
the critical point, corresponding to each distribution. Fi- 
nally, the conclusions are summarized in section V. 

II. THE HAMILTONIAN AND THE 
MIGDAL-KADANOFF 
RENORMALIZATION-GROUP PROCEDURE 

Let us consider the nearest-neighbor short-range Ising 
spin-glass model on a general p-connection diamond hi- 
erarchical lattice. The hierarchical lattices are connected 
graphs, recursively constructed by replacing all bonds at 
each generation by a basic unit, the starting point be- 
ing the basic unit itself (first generation). The DHL ba- 
sic unit is designed by two root-sites coupled through p 
connections in parallel, each of them composed by two 
bonds in series via one internal site (scaling factor 2), as 
schematically shown in Fig. |l|. The graph-fractal dimen- 
sions of such lattices are given by d = 1 + In p/ In 2. 




FIG. 1. General diamond hierarchical lattice basic unit, 
with p parallel paths, each one composed by two bonds and 
one internal site (scaling factor 2). 

The Hamiltonian of the present model is given by 

-m= K N (i,j)o- i0 - h (i) 

<i,j> 

where Kpj(i,j) = — /3Jjv(«,j) is the reduced exchange 
coupling constant between nearest neighbor spins of a 
A-generation DHL, Jjy(i,j) being the N th -step Migdal- 
Kadanoff renormalized coupling constant obtained from 
the original distribution. It is worth to mention that the 
real-space Migdal-Kadanoff renormalization-group trans- 
formation on cZ-dimensional Bravais lattices is known to 
be equivalent to the exact solution on the DHL with a d- 
graph fractal dimension [ fL4| . The MKRG transformation 
of a p-connected DHL with general non-uniform coupling 
constants yields an effective coupling given by pT| , 
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«n U*, M ) - 2 2j I cosh [if M ) - tffo ; ^ J ' 

(2) 

which can be rewritten in terms of more convenient vari- 
ables as 



t^i = tanh 



i=i 



(3) 



where t^ = tanhif^^ is called thermal transmissivity. 
In Eq. (g) ai labels the spins variables of the p internal 
sites of a basic unit while /x and /x' denote the root ones 
(see Fig. |). 

The critical temperature for the spin-glass model can 
be numerically obtained either from eqs.(|2|) or (|3|) by 
monitoring the width of the iterated distribution 
For the case of Eq. (^), in the zeroth step a pool of 
M random initial transmissivitics ({ij} = tanh(/3{ Ji}); 
i = 1,2,..., M) is generated, with the coupling constants 
{Ji} following a given probability distribution. In the 
first step Eq. @ is iterated M times by choosing at ran- 
dom 2p initial transmissivities; the new M-valued pool 
represents the renormalized thermal transmissivity dis- 
tribution. This process may be repeated and renormal- 
ized distributions can be numerically followed by comput- 
ing its moments at each step jllj and a flow diagram may 
be constructed, e.g., in the transmissivity versus variance 
plane @. Since the spin glass fixed-point distribution is 
not analytically known, we considered four initial sym- 
metric distributions of interest, namely the Gaussian, Bi- 
modal, Exponential and Uniform ones, defined below, 



P(Jij) 



s/2tt 



exp(- 



h Jfj) (Gaussian), 



(4a) 



P(Ji,j) = 3 [8(Jij ~ 1) + S(Jij + 1)] (Bimodal), (4b) 
P(Ji,j) = 75exp(- v / 2|Jij|) (Exponential), (4c) 



and 







otherwise 



(Uniform) . 
(4d) 



To each initial distribution a pseudo-critical tempera- 
ture is associated, for which the flow will converge to the 
critical point characterized by the " fixed-point" distribu- 
tion which is numerically known [ pl| . For temperatures 
very close but below (above) this pseudo-critical temper- 
ature the flow will at first, approach the "fixed point" 
and then turn on to the spin-glass (paramagnetic) fixed 
point characterized by infinite (zero) variance and zero 
mean. In Table 1, we present our numerical estimates 
for these pseudo-critical temperatures obtained from the 
initial distributions given by eqs.(Q) in the case d = 3 
(p = 4), compared to other values reported in the litera- 
ture. 



III. THE LOCAL MAGNETIZATION AND THE 
LOCAL EA ORDER PARAMETER 

A. The generalized method for spin glasses 

The aim of this method is to establish a recursion re- 
lation between the value of the local magnetization of 
the internal site belonging to a connection of a given ba- 
sic unit of the DHL with the local magnetization of its 
root sites. If this is achieved, one can consider a finite 
DHL with N generations with Ising spin variables, de- 
scribed by the Hamiltonian given in Eq. (|l|), with nearest- 
neighbor random exchange coupling constants chosen 
from a given initial distribution, and renormalize it N— 1 
times storing at each step all renormalized coupling con- 
stants. Then, taking arbitrarily initial magnetizations 
(corresponding to the spin-glass boundary conditions) for 
the root sites of the first generation, we can successively 
calculate the local magnetizations of each site of the N- 
generation DHL, at a given temperature and for a chosen 
initial coupling-constant distribution. 

To obtain this recursion relation, let us consider a N- 
generation DHL and look at an arbitrary basic unit in- 
troduced at the TV^-generation as shown in Fig. [| This 
basic unit is connected to the lattice by its root sites 
(/x, //'). Therefore, the partition function of the whole 
lattice can be written as 

Z = Tr ({CTiW) exp {-/?#'} 

= Tr^^^vi-pHiio-i},^^)} (5) 
x exp{— (3 [hftfi + h/j/fj + K' /x/x']} 

where {oi}, i = l,2...p, denote the internal spins 
within each connection, /x, p! are the root spins of the 
basic unit, h^, h^i and K' are, respectively, the effec- 
tive fields and the effective coupling acting upon the ba- 
sic unit root spins induced by the rest of the lattice, 
i? ({<7i}, /i, //) is the internal Hamiltonian of the basic 
unit given by 



p 



H({o- i }, f i,n') = ^(K it i- 



(6) 



where Ki, K' { are the corresponding random coupling con- 
stants between the oVs and the /x, /x' spins, respectively. 
The local magnetizations of all sites within a given basic 
unit can be easily evaluated for the effective model 

TABLE I. Spin-glass pseudo-critical temperatures for dif- 
ferent probabilities distributions within the Migdal-KadanofF 
renormalization-group procedure, in d — 3 spin glass. 





Gaussian 


Bimodal 


Exponential 


Uniform 


This 
Ref. [ 
Ref. [ 
Ref. [ 
Ref. [ 
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H 

20 


rk 

1 
1 
1 
1 


0.88 
0.88±0.02 
0.85 
0.89 
0.83 


1.15 
1.05±0.02 
1.2 

1.15 


0.75 
0.7 
0.71 


0.96 
1.00±0.02 
1.0 

0.94 
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defined by the Hamiltonian H' by tracing over all spins 
variables. However, since our main concern is to establish 
a recursion relation between the internal-site magnetiza- 
tion and these of the roots sites of a certain basic unit, the 
procedure can be further simplified by focusing our atten- 
tion to a single connection and including the effects of the 
other connections onto the effective fields and coupling, 
following the ideas of the decoration transformation for- 
malism p2| . In this case, our system is over- reduced to a 
single connection with an internal site under the action 
of the effective fields and effective coupling induced by 
the remaining lattice, as schematically shown in Fig. |[ 




FIG. 2. Schematic representation of the equivalent system 
constructed by retaining a given basic unit of the final gener- 
ation, with coupling constants Ki (i — 1, 2p). fj,, \i! are the 
root spins of the basic unit, h^, h„i and K' are, respectively, 
the effective fields and the effective coupling generated by the 
whole lattice. 

For this over-reduced system, the Hamiltonian in de- 
fined by 

H" (a, (j,, fjf) = (K-lh + K 2 fJ.')a + h^i + Vm' + K ' MA*' • 

(7) 

The magnetizations at each site are straightforwardly cal- 
culated, giving 

< n > = Z~ x Tr \jj,exp(-0H")} 



1 + tTpTtf + t x t 2 {t + r M r M /) ' 
<fjf>= Z~ l Tr [fjf exp(-/3H")} 

_ (y + 1^) + ti^fa + fry) 

1 + tr M T^ + ht 2 (t + t m t m >) ' 

< a > = Z~ l Tr [aexp(-f3H")] 

_ hXjp + fry) + f 2 (y + trjt) 
1 + fr^v + ht 2 (t + yy) ' 



(8) 



(9) 



(10) 



where t = tanh(K ), tj = tanh(iQ) (i = 1,2) and r Q = 
tanh(/i Q ) for a — fJ,,fi'. Now, from Eqs. (g) and (||) 
we can write the unknown variables (i, and r M ') as a 
function of < a > and < fjf > and Z, and substitute 
them in Eq. (jlO|) to end up with the recursive equation 



< a >= 



i 



f 2 f 2 



< fJ, > 



\-t\t\ 



< /i > 



(ii) 



We emphasize that if the sites of a hierarchical lattice 
are properly addressed, Eq. establishes a recursive 
equation between the local magnetization of the sites be- 
longing to the last generation and the ones of previous 
generations. Moreover, its coefficients do not depend on 
the unknown fields and couplings, but only upon the cou- 
pling constants belonging to the chosen connection. This 
result is the main achievement of this method. 



B. Numerical Procedure 



To analyze the structure of the local EA order param- 
eter of our model, we should calculate < oi > 2 for all 
sites and average them over many samples, yielding to 



-EA 



= < a, 



> 2 1 



(12) 



where [...] c stands for the configurational average taken 
over many independent initial distributions of couplings. 
To consider larger lattices we have to go further in the 
renormalization steps. Since the number of sites and 
bonds increases like (2p) N , the amount of computer mem- 
ory required to store the magnetizations and the coupling 
constants during the intermediate steps will increase with 
such rate. In order to maximize the number of renormal- 
ization steps, we look at the magnetization structure of 
a subset of representative sites of the lattice. These sites 
are the 2 N ones belonging to any shortest path connect- 
ing the roots sites. The magnetization (and/or the EA 
local order parameter) structure of this subset can be 
viewed as a representative profile of the whole lattice. 
Since they are stochastically equivalent, we argue that 
by averaging over many profiles, we should obtain the 
correct scenario for the local EA order parameter of the 
considered model. 

To calculate the profile of the local magnetization of a 
N generation hierarchical lattice we make use of Eq. ([ll]). 
To display the profile, we have to label the sites of a given 
path, assigning the values of its local magnetizations of 
the support set, due to the graph topological nature of 
the hierarchical lattices. To proceed, we choose the set of 
site labels by (s, I) belonging to the interval [0, 1], defined 
by s x 2~ l , s = 1,3,5, {2 l - 1), and I labeling the 
generation (I = 1,2,..., TV). For this choice the recursive 
equation can be written as 



< a > s ,i= A Suh < fi > Sl ,i 1 +K S ,U < M' >s 



(13) 
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where si = |(s ± 1), l± = I — 1, Sj — his^fl), lj = I — j, 
with j — 2, 3, I. The coefficients of Eq. ( |13| ) are given 
by 



A. 



s i.'i 1 _ +2 +2 

- 1 - 1 'S,Sl''S,Sj 



A 



t.^l-t*^) 

1 _ ^S,Si^S,Sj 



(14a) 



(14b) 



where t„ 



= tanh[AT;(,s, sj)], Ki(s, Sj) being the coupling 
constant between the spins at the positions s.2~~' and the 
one at Sj.2~ lj . 

To calculate an EA order parameter profile we must 
first generate the coupling constants for each level. Due 
to the disordered nature of the profile we make use of an 
equivalent stochastic procedure in order to save computer 
memory at intermediate steps of the calculation. For a 
fixed value of the temperature we create an initial distri- 
bution for the thermal transmissivities, represented by a 
pool of M random numbers (M ^ 10 • 2 N ). At the Nth 
level we choose randomly, from the initial distribution, 
a set of 2 N couplings, which are stored to be used later 
in the calculation of the site magnetizations. At the next 
level [(N — l)th level], we obtain a renormalized distri- 
bution (new M random numbers), generated according 
to the renormalization Eq. (||); from this distribution, 
we pick randomly 2 Ar_1 couplings, which are also stored. 
This process is carried for N — 1 times, such that at the 
last level, only two couplings are stored. Now, we make 
use of Eq. (p3|), fixing the initial values for the magneti- 
zation of the roots (zeroth generation) and calculate the 
local magnetization of each level, using for the coupling 
constants the values previously stored. 



EA 



(a) 
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3x10 



Position 



4x10 



FIG. 3. Portion of the profile of local EA order parameter 
for one sample and a lattice with N = 16 generations; the 
sites chosen belong to a sub-set with positions in the range 
i = 3.0 x 10 4 , 4.0 x 10 4 , selected from the set i = 1, 2, 2^, 
corresponding to a given shortest path connecting the two 
roots sites, (a) T = T C . (b) T = 0.9T C . 



IV. THE MULTIFRACTAL PROPERTIES 



C. The EA order parameter profiles 

For each of the initial distributions defined in Eq. 
(||), we generated profiles at the corresponding critical 
temperature T c (see table |), as well as T± = 0.9T C , 
T 2 = 0.8T C and T 3 = 0.7T C . This was done for lattices 
with d = 3 and N — 8 up to 16 hierarchies. In Fig. ||| we 
display these profiles for the Gaussian distribution at the 
temperatures T c and T± (N — 16) for just one sample, 
whereas in Fig. |], the same is done for 200 samples. It 
is clear from these figures that the disordered structure 
of the local EA order parameter of one sample increases 
as we go further in the condensed phase. However, when 
the configurational average is taken, the profile presents 
uniformities reminiscent of the graph lattice symmetry, 
similar to what happens for the pure model [ jl5| . It is also 
evident from Fig. [|, the increasing of q EA = 2~ N ^ i qf A 
, the mean value per site, as we decrease the tempera- 
ture. For all other distributions listed above (Bimodal, 
Exponential and Uniform), the same qualitative behavior 
was observed [E3|l . 



The high degree of discontinuity shown in the profiles 
suggests us the multifractal analysis as a tool to investi- 
gate the singularities of the measure constructed from the 
EA local order parameter, following the same approach 
used for pure models jF5] [17| . To obtain the multifractal 
spectra (F (a )- function), we first define a measure by the 
normalized local EA order parameter 



c, 



EA 



^EA 



Ei<2 



EA ' 



(15) 



and construct a parametrized family of normalized mea- 
sures defined by 



EA 



(C 



EA\q 



(16) 
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FIG. 4. Portion of the profile of the average local EA order 
parameter (200 samples) for a lattice with N = 16 genera- 
tions; the sites chosen are the same as in Fig. |^. (a) T = T C . 
(b) T = 0.9T C . 





a 



The FfaJ-function is now obtained following the 
method due to Chhabra and Jensen (24) , where the spec- 
trum is obtained by varying the parameter q and calcu- 
lating 

In Fig. |5|, we display the corresponding i*Ya,)-functions 
for the profiles of one sample obtained from the initial 
distributions listed in Eq. (0), for temperatures at and 
below T c . We notice light variations of the spectra by 
changing the temperature for the Gaussian, Bimodal and 
Uniform cases, whereas more pronounced changes are ob- 
served for the Exponential case. In Fig. [| the corre- 
sponding Ffaj-functions averaged over 200 samples are 
shown for the Gaussian and Bimodal distributions at T c 
and Ti = 0.9T C . Minor changes are observed when we 
compare with Fig. g. 



FIG. 5. F(a)-functions of the local EA order parameter 
profile with one sample, for decreasing temperatures and dif- 
ferent initial distributions of coupling constants, (a) Gaus- 
sian, (b) Bimodal, (c) Exponential and (d) Uniform distri- 
butions (O T = T c , □ Tl = 0.9T C , O T 2 = 0.8T C and A 
T 3 = 0.7T C ). 

In Fig. ^, we show in the same plot the i^aj-functions 
at and below T c for the four considered distributions. It 
is worth to call the reader's attention to the universal 
character of the F (a)-ixmction within small deviations. 
Here we remark that although we have used the renor- 
malized distributions of couplings at each step of the cal- 
culation, the influence of the initial distribution should 
be relevant, since in the thermodynamic limit (N — * oo), 
half of number of sites in the profile belong to genera- 
tion N and their magnetizations are calculated with the 
coupling constants introduced by the initial distribution. 
For the whole lattice this influence should be even more 
relevant; in this case | of the total number of sites belong 
to generation N. 



F(a) 





a 



a 



FIG. 6. F(a)-functions of the average local EA order pa- 
rameter profile (200 samples), at and below the critical tem- 
perature, for the (a) Gaussian and (b) Bimodal distributions 
of coupling constants (Q T = T c , □ Ti = 0.9T C ). 



F(0C) 





(b) 



the temperature, for T > T c . On the other hand, the 
a m i n exponent which governs the set of higher measures 
should remain finite to describe the singularities of the 
measures belonging to the sites "close" to the root sites 
(or surface sites). Those sites are the ones whose magne- 
tizations were calculated with at least one of the values 
imposed as initial boundary conditions. The intermedi- 
ate points of the spectra should be spurious points, since 
the present algorithm Q used to calculate the F(a)- 
function is based on the method of moments and tend 
to produce a top envelop of the actual spectrum |25|] . In 
the thermodynamic limit (N — > oo), we expect to have 
no spectrum except a point (0, 0) corresponding to the 
non-vanishing values introduced by the imposed bound- 
ary conditions. 
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FIG. 8. Upper (O awx) and lower (□ a m in) bounds of 
the F(a) spectrum as a function of the temperature for dif- 
ferent distributions of coupling constants, (a) Gaussian, (b) 
Bimodal, (c) Exponential and (d) Uniform distributions. 



FIG. 7. Comparison between the _F(a)-functions of the lo- 
cal EA order parameter profiles at and below the critical tem- 
perature for the O Gaussian, □ Bimodal, O Exponential and 
A Uniform initial distributions of coupling constants, (a) 
T = T C and (b) Ti = 0.9T C . 

In order to investigate the behavior of the F(a) spec- 
trum around T c , we plot in Fig. || the dependence of 
its upper and lower bounds as a function of the tempera- 
ture, for all the considered initial distributions. We notice 
that below T c a constant value is observed for the upper 
bound (a max ), while small deviations occur for the lower 
bound (a m i n ) close to T c . Nevertheless, for temperatures 
slightly above T c an abrupt increase is observed for the 
upper bound while a small decrease occurs for the lower 
bound, signalizing the SG transition. At T > T c the mag- 
netization at the majority of sites vanishes, being elimi- 
nated from the calculation of the spectrum. Therefore, as 
we are dealing with a finite lattice one should expect a fi- 
nite but higher a max , which is the exponent governing the 
singularities of the set of smallest measures still present. 
This is evidenced by the rapid increasing of a max with 



V. CONCLUSIONS 

We generalized the exact recursion method developed 
by Morgado et al. |15j and applied it to investigate the 
structure of singularities of the local Edwards- Anderson 
order parameter of the short-range Ising spin-glass model 
on diamond hierarchical lattices. Within this procedure, 
the distribution of coupling constants is renormalized by 
the Migdal-Kadanoff renormalization equation, which is 
exact in such classes of lattices. The EA order param- 
eter profiles (2^ sites) for lattices up to N = 16 hier- 
archies were calculated by considering four types of ini- 
tial distributions (Gaussian, Bimodal, Exponential and 
Uniform) and for temperatures around the critical point. 
For T <T C the profiles of the spin-glass condensed phase 
show a high degree of disorder and singularities increas- 
ing for lower temperatures. The multifractal analysis was 
applied to these profiles, revealing a large spectrum of ex- 
ponents for the singularities of the measure defined by the 
normalized local EA order parameter. For each consid- 
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ered distribution of couplings, the F(a)-function of the 
profiles show slight variations with the temperature in 
the studied range of 0.7T C ~ T c , with larger deviations 
occurring when the temperature gets closer to T c . These 
profiles reveal a high degree of local disorder of the spin- 
glass condensed phase, a scenario not observed in the 
pure case [ fulfill for the same classes of lattices. More- 
over, small deviations are observed when the multifractal 
spectrum obtained from distinct initial distributions of 
coupling constants are compared at the same tempera- 
ture. This suggests an universal multifractal behavior of 
the present model for distinct initial distributions of cou- 
plings, taking into account that | of the total number of 
sites belong to the last generation, whose local magneti- 
zations are calculated with the coupling constants intro- 
duced by the not-yet-renormalized (initial) distribution. 
We have also studied the temperature dependence of the 
range of the multifractal spectrum close to the critical 
point. For T < T c the range of the a- Holder exponent 
remains quite constant. However, when the temperature 
is higher than T c , one observes an abrupt change in the 
multifractal spectrum signalizing the transition. 

Contrary to the pure case, where a non-trivial mul- 
tifractal behavior is observed only at the critical tem- 
perature p^pll , the persistence of the F(a)-function 
throughout the spin-glass phase indicates the highly non- 
trivial character of such phase. Although we are not able 
to associate the persistence of multifractality with any 
prediction from the available theories to describe short- 
range spin glasses, the present work evidences the con- 
trast between the spin-glass and ferromagnetic states and 
the critical nature of the spin-glass condensed phase. 
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